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Abstract We develop a nonlinear, three-dimensional phase field model for crys¬ 
tal plasticity which accounts for the infinite and discrete symmetry group G of 
the underlying periodic lattice. This generates a complex energy landscape with 
countably-many G-related wells in strain space, whereon the material evolves by 
energy minimization under the loading through spontaneous slip processes induc¬ 
ing the creation and motion of dislocations without the need of auxiliary hypothe¬ 
ses. Multiple slips may be activated simultaneously, in domains separated by a 
priori unknown free boundaries. The wells visited by the strain at each position 
and time, are tracked by the evolution of a G-valued discrete plastic map, whose 
non-compatible discontinuities identify lattice dislocations. The main effects in the 
plasticity of crystalline materials at microscopic scales emerge in this framework, 
including the long-range elastic fields of possibly interacting dislocations, lattice 
friction, hardening, band-like vs. complex spatial distributions of dislocations. The 
main results concern the scale-free intermittency of the flow, with power-law ex¬ 
ponents for the slip avalanche statistics which are significantly affected by the 
symmetry and the compatibility properties of the activated fundamental shears. 

Keywords crystal plasticity; intermittency; dislocations; finite strain; phase field 
modeling 
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1 Introduction 

The high-resolution investigation of plastic flow in crystalline materials at mi¬ 
croscopic scales shows, in contrast to its smooth macroscopic behavior, a bursty 
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dynamics related to the nucleation and motion of lattice dislocations, with power- 
law statistics for the slip avalanches [2ll[9l[571|63l|26]. Theoretical work on these 
phenomena has been developed by using discrete dislocation dynamics [371II9], 
pinning-depinning |39| , or automaton models |491l5n| , renormalization group meth¬ 
ods m, and continuum descriptions based on Nye’s dislocation density tensor 
|25| . Phase field methods have also provided fruitful tools for the analysis of many 
aspects of crystal plasticity [58l[5^ll71l30ll31ll5Tll60ll35l[3^[33ll36ll3l| . 


As defect nucleation, annihilation, and movement originate from concentrated 
slip processes which accommodate large strains while locally striving to preserve 
the original lattice structure, here we propose a new nonlinear three-dimensional 
(3D) phase field model accounting for crystal symmetry. This is described by the 
global symmetry group G of crystals, conjugated to the group GL(3,Z) of uni- 
modular 3x3 integral matrices preserving 3D Bravais lattices [23l[6l[Ti|l3|l43] . The 
group G dictates the location and distribution of countably-many relaxed states in 
the space of strain tensors (see also a related approach in [18]). In this framework 
the plastic distortion map is therefore discrete as it is G-valued, parameter¬ 
izing the energy well locally visited by the material’s deformation gradient F at 
any x and t. The discrete jumps in the plastic map F^ are related to the presence 
of lattice dislocations in the material, whose long range elastic fields agree with 
the classical linear elastic predictions, emerge as non-compatible discontinuities of 


Near any F^ £ G the amount of shear along each principal slip direction of the 
crystal constitutes a scalar phase field component. These fields select the deforma¬ 
tion path through energy minimization, and the coupling with the discrete plastic 
distortion allows the material system to explore the complex energy landscape 
induced by crystal symmetry. The evolving phase fields and body deformation 
lead in turn to a crystallographically-informed dislocation dynamics, driven by 
energy relaxation, which automatically accounts for the possible presence of si¬ 
multaneously activated, interacting multiple slip modes, without the addition of 
supplementary hypotheses. 

The model can resolve isolated defects in the crystal as well as configurations 
possibly involving complex defect topologies. By simulating the behavior of the 
material under cyclic mechanical loading, we show that this framework describes 
correctly various features of the plastic flow in crystalline materials, and sheds 
new light on the properties of the associated scale-free intermittent behavior at 
microscopic scales. 

The paper is organized as follows. In Sect. we describe the phase field model, 
including the free-energy functional driving the gradient-flow dynamics, together 
with the new discrete flow rule governing the evolution of the plastic distortion F^. 
In Sect. 1^ we evidence how F^ characterizes dislocations within the material, and 
we test the model by analyzing numerically the mechanical response of a single 
crystal subject to different types of cyclic loadings. We also study some of the 
statistical properties of the bursty dynamics emerging in these simulations. We 
draw the main conclusions in Sect. |4| 
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2 Model 

Let an initially homogeneous stress-free single crystal occupy a reference domain 
fl in 3D space. We denote by y the displacement field, with F = Vy the de¬ 
formation gradient, J = detF, and, respectively C = F^F and E = ^(C — I), 
the Cauchy-Green and strain tensors. Let = e^jU^ (summation understood, 
with j, h = 1,2, 3) be a basis for the reference crystal lattice, where Uj are pre-set 
mutually-orthogonal unit vectors. The global symmetry group is then given by 
G = E~^GL{3,Z)E, with E = (e^j), see [l3] . 

Further information regarding the crystallographic features of plastic flow in 
the material is given by the p principal slips, defined by the shears Sq, = I ± ba (g) 
tia G G, mapping the crystal to itself (a = 1,... ,p). The Schmid tensors bo, ® Hq 
are obtained from mutually orthogonal vectors ba,na!, giving respectively the a- 
th slip direction and slip plane normal. These are dictated by the crystallographic 
and mechanical properties of the material. 

We assume the standard Kroner-Lee decomposition of the deformation gradient 
F = FeFp in terms of an elastic and a plastic component, where Fp parameterizes 
the relaxed states of the material Ezma. As the latter are given here by the tensors 
in G, we are led to the natural assumption that Fp coincides with the discrete field 
F^ G G, defined as the bottom of the G-related energy well to which the gradient 
F belongs at any given place x and time t. We further assume the elastic distortion 
Fe to be given by a component Fp,, accounting for the activated principal lattice 
slips, multiplied by a possible further residual distortion Fr. Explicitly, we set 

F = FeFp = F.F^Fp^, (1) 

with 

p 

Fp = FpGG and F^ = I-|-^ 0 (pq) b^ (g) n^, (2) 

a=l 

where the variables pa G [—1,1] measure the amount of shear along each Sq. 
As in [35], the role of the function = 3s|s| — 2s^ in ([^ 2 , which is odd and 
monotonically increasing in [—1,1], with stationary points in 0 and ±1, is to locally 
extend the stability range of all the G-related strain equilibria considered below. 

Due to assumptions 0-§, here the variables pa act as phase fields interacting 
with the deformation field y. This is further implemented by the assumptions we 
make for the anisotropic free energy density tj) of the material, which is taken 
to have the form ■!/) = Tpp + Tpr, where the phase field related term tpp has the 
piecewise-polynomial expression 

tpp{p) = Ai ^Pa(l - \pa\)^ +A 2 pip}, (3) 

a 

with p = (pi,..., pp), while tpr gives an isotropic penalty, with elastic modulus K, 
to the residual distortion Fr: 

l/>r = ifflFr^Fr-ip. (4) 

For Afi > 0, 'ipp is minimized when at most one among the pa is equal to the 
transition value ±1, and the other vanish [53]. The potential ipp enforces the re¬ 
quired symmetry under permutations of the pa, and also implies the symmetry 
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with respect to the backward and forward motion through the principal shears. 
The ^ 2 -term favors the activation of one single slip at a time. The ratio A^jAx 
determines the width of the optimal shearing paths close to the coordinate axes 
whereon the system moves from one minimizer to another, as the greater An the 
smaller the remaining phase fields may be when one of them departs significantly 
from zero. Finally, the ratio AxjK controls the height of the energy barriers be¬ 
tween slip-related minimizers (and therefore the plastic yield stress). The adopted 
minimal expressions for ipr and xlxp, model two dominant effects occurring in plastic 
shear, i.e. the presence of (tunable) energy barriers, and the preference for sin¬ 
gle slips. A more precise characterization of ip would allow to mimic the detailed 
structure of the energy landscape of crystals within each well, as with anisotropic 
7 -surfaces given for instance by ab-initio models |38II121I20] . However, we keep the 
constitutive choices in the present framework as simple as possibl^ 

The evolution of the pair (y, fi) is determined by a gradient-flow equation of the 
Ginzburg-Landau type, coupled with the standard Cauchy equation for mechanical 
equilibrium 


dfi _ 6ip 
dt 5fi ’ 


diver = 0, 


(5) 


where er = is the Cauchy stress. As with to limit the num¬ 

ber of constitutive parameters in the model we consider in an isotropic mobility 
coefficient L multiplying the variational derivative of ip with respect to /x. Equa¬ 
tions 0 are solved subject to initial conditions for /x, and to either displacement 
or traction boundary conditions for y on dO. 

Along the evolution of /x, the map must be discretely updated at any t > to 
when one of the phase fields, say /xa, reaches the value ±1, as this indicates the 
local attainment by F of a G-neighboring relaxed state under the driving: 


Fp^ ^ (I±ba(g)n„)Fp^, (6) 

the corresponding pa being reset to zero. Rule 0 provides a tensorial general¬ 
ization of the scalar updating considered in previous simplified models of crystal 
mechanics [42114911^ . The remaining phase fields, other than pa, are updated in 
such a way to guarantee the continuity of the combination F^^F^ in 0 across jump 
instants. With F continuous by hypothesis, the above means that Fr, ip, ipfi, ipr and 
a are also always continuous, despite the discontinuities inherent in the evolution 
of the discrete map F® in 0 and in the corresponding updated fia- An exception 
to this general behavior occurs when, during the evolution of /x, the resulting F^^ 
in 02 is not isochoric. For instance, this happens when for some a, pi, two shears 
Sa and are active in ([ 2|2 at some jump instant, with (bp • n^)(b^ • na) i=- 0 and 
non-zero (pit^a) and (p{ii^. In this case, if pa reaches a transition value ±1, then 
fip is set equal to zero to recover the natural isochoricity assumption on F^^. The 


^ As for instance in m, we do not include here any interface-penalizing gradient term in the 
free-energy functional because in all our simulations we keep the typical distance H between 
interpolation nodes always greater than or equal to the minimal physically relevant length h 
in the continuum, which in the present context is taken of the order the lattice cell spacing. To 
model plastic behavior of crystals at scales H larger than h, besides including a gradient term 
in the phase field potential, the Schmidt tensors in 0 must also be suitably re-scaled (see for 
instance ESI). The corresponding scaling imposed on the potential determines whether in 
the limit H ^ h the yield and Peierls stresses vanish or stay finite. In what follows we only 
consider plasticity at microscales with ~ /i. 
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resulting discontinuity of Fr may be kept under control, as it becomes smaller the 
larger we choose A 2 in t/i^. The cross term proportional to ^ 2 , indeed, penalizes 
non-zero values for any /x components when any of them approaches ±1. 

The discrete updating rule ^ for has here the role that is played in the 
classical continuum theories of plasticity mm by the standard differential flow 
rule for Fp, which is used also in phase held models for martensitic transformations 
and plasticity in crystals |35U321I331I34| . However, represents a signihcant novelty 
of the present approach, as it enforces the condition that the material locally 
reaches a stress-free equilibrium through only if F belongs to G modulo 

an orthogonal transformation. This is because, in an unperturbed system F is a 
solution ofra only when both the phase held vector /x minimizes V'p (and therefore, 
from Q-([2|7FpF^ £ G), and Fr is orthogonal, as required. In previous phase held 
implementations of crystal plasticity (either in the linear |581l591ll71l31U51ll36| or in 
the hnite-strain |351l32l[33ll34| regimes), the relaxed states obtained by using the 
classical differential how rule result to be lattice invariant maps only when a single 
slip is active because only in this case integer values of the phase helds correspond 
to tensors in G. By contrast, in the general case the values of Fp obtained from 
the classical how rule bear no relation to the underlying lattice structure of the 
crystal expressed by G. 


3 Results 

Dislocations. The discrete held F^ has a crucial role in locating dislocations, 
which are associated with its non-compatible spatial jumps. We recall that in clas¬ 
sical continuum theories m the plastic distortion map Fp describes a dislocation 
whenever it lacks compatibility, i.e. when it is not possible to write it as the gradi¬ 
ent of a vector held. In the present setting based on the discrete plastic distortion 
F®, dislocations may occur in regions where the map F® is not compatible (ana¬ 
lytically, where curlF^ 0, though this should be intended in the distributional 
sense as F^ is a piece-wise constant held). As a typical instance, we consider a 
partition of the reference domain 17 into subsets on which F^ is constant. The 
induced jump discontinuities of F^ are compatible if they satisfy the Hadamard 
conditions |[F®]] = b (g) zx, with u the unit normal to the jump surface [33]. This in 
particular implies that any compatible discontinuities must necessarily occur along 
(portions of) planar surfaces. When two such parallel surfaces, together with a fur¬ 
ther surface connecting them (which is necessarily non-compatible), are considered 
on a small neighborhood of a line, this latter is identihed as a dislocation line in 
the lattice, so that compatibility failures - including the presence of non-straight 
portions of jump surfaces - may induce surface dislocations [TIITOUT] . 

An explicit example illustrates how dislocations are described by the non¬ 
compatible discontinuities of F^. We consider a simple-cubic crystal in 17, with 
mutually-orthogonal lattice vectors ej = auj, j = 1,2,3, where a is the lattice 
parameter, and Uj are unit vectors along the 4-fold axes. In this lattice there are 
p = 6 fundamental slips Sq , which involve all the mutually-orthogonal couples of 
vectors u_,-. In this case the Sq, generate multiplicatively the entire symmetry group 
G, which for this Bravais type coincides with GL(3,Z). 

As in Fig. we select 17' C 17 with boundary coinciding with the jump surface 
for F®, composed by two nearby planes orthogonal to 62 , with distance h (see 
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Fig. 1 (Color online) Left: Computed stress field near a dislocation. Right: Comparison of 
the results obtained for as a function of the distance from the center of the dislocation, 
for X = 0: present model (solid line) and linear-elastic solution (dashed line). 


footnote and connected by a (necessarily non-compatible) portion of surface 
close to a (dislocation) line yais- By considering now a closed curve 7 in 17 inter¬ 
secting the jump surface in two points Pi and P 2 , one on each parallel plane, with 
P 2 = Pi -|- he 2 , we compute dS = |[F®]](PiP 2 ) = hh, with t the unit tan¬ 

gent along 7 . Therefore, in this case hh is the classical Burgers vector associated 
with the dislocation EZl. The standard edge and screw dislocations are retrieved 
as non-compatible discontinuities of F® for which b is respectively orthogonal or 
parallel to the line 7 dis. 

To compare the stress held computed for a dislocation with the classical results 
derived from linear elasticity, we set F^ to be respectively equal to I and I-|-ui®U 2 , 
in and 17'. An edge dislocation results in 7 dis, with Burgers vector hui. The 

(rxx-component of the corresponding stress held is plotted in the left panel of Fig.[^ 
computed with the Comsol - Finite Elements Method software, used hereafter in all 
numerical simulations. In the right, this solution is compared with the analytical 
prediction [28], with shear modulus and Poisson ratio dictated by the potential 
We see that the numerical simulation matches the analytical prediction to 
distances of the order of the defect width h. 

Starting from a local minimizer identifying a lattice dislocation as in Fig. 
the Peierls stress ap is given by the external load at which the defect hrst becomes 
unstable, and here can be computed by considering a monoparametric family of 
equilibria of the system where all the phase fields are quenched except for the 
one associated with the principal slip active in 17'. The value of ap depends on 
all the energy constants, and by setting (here and in all the ensuing numerical 
simulations) the elastic modulus to K = 200 GPa, with AijK = 1, A 2 /K = 5, and 
h = 1 nm, we obtain crp ~ 1.6 GPa, in line with values computed in other theories 
of crystal plasticity IHIIHI. 

Plastic flow under cyclic loading. We now study the response of a single 
crystal to a possibly cyclic loading to show explicitly how the creation and evolu¬ 
tion of dislocations spontaneously progresses during plastic flow in this model. We 
consider a system initially in local equilibrium, with vanishing phase field compo¬ 
nents, and Fp = F^. The well-tracking map F^ is initially set to coincide with the 
identity, except for the presence of possible quenched defects in the body, modeled 
by randomly distributed small regions where the phase field components and the 
map Fp ^ I are not allowed to evolve. We impose hve different loadings I-V, on a 
simple-cubic crystal as above, occupying a box domain 17 with free lateral sides. 
The bottom of 17 has unit normal m, and the lattice basis ej is along the 4-fold 
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Fig. 2 (Color online) Intermittent plastic flow under cyclic shearing, in a box domain Q 
with two edges parallel to a and m, containing few randomly-positioned quenched defects. 
Relevant components of the stress and strain tensors are monitored. Left: Serrated stress-strain 
hysteretic loop for a typical cycle of loading-unloading with a single active slip (simulation I), 
and no hardening. Inset: decrease of the stress peak on the first loading, for an increasing 
number of identical defects randomly distributed in the sample. Right: Stress-strain curves 
in other simulations. Loading II (continuous blue and red dotted lines) mostly involves two 
compatible activated slips. No hardening is observed as for case II. The results of simulation 
III (not shown) also involve mostly compatible slips and produce no hardening, analogously to 
cases I and II. Simulation IV (green), and simulation V (purple) mostly involve non-compatible 
slips and show a clear hardening effect. 


axes. In simulation I the basis ej is aligned with all the sides of 17, while in simu¬ 
lation II the crystal is rotated around the bottom normal m by about 45 degrees. 
The boundary conditions fix the bottom of 17 and impose a top-face deformation 
consistent with an overall shear I + A a (g) m, where the vector a is orthogonal to m 
and parallel to one of the free box faces. These boundary conditions are designed 
to activate mostly one or two (compatible) slip systems during plastification. In 
simulation III the same boundary conditions are imposed but the lattice basis in 
17 is rotated by about 30 degrees around the vector a, so as to possibly activate 
more than two compatible slip system in the lattice. The more general loadings IV 
and V are then imposed to activate mostly non-compatible slips. Specifically, in 
case IV the isochoric deformation on the top and bottom faces of 17 is consistent 
with I -|- A(ei (g) 63 -|- 63 (g) 62 ) , while in simulation V an arbitrary rotation is con¬ 
sidered of crystal axes ej in the box, with the same imposed boundary conditions 
as in cases I, II, III. The parameter A is driven quasi-statically, so that the system 
evolves by local minimization through inhomogeneous equilibrium states, obtained 
by solving the equations in coupled with (§. The resulting deformation fields 
are expected to activate a single principal slip in simulation I, and (mostly) two 
kinematically-compatible principal slips in simulation II (thus allowing in princi¬ 
ple for defect-free interfaces). In contrast, in simulations III, IV and V, several 
slips are activated, possibly compatible or non-compatible among each other. In 
all cases the map F® tracks where and how such slips progress and eventually 
interact in the material, with the finite and multiplicative updating rule for 
F® guaranteeing the relaxed configurations to locally be driven to those dictated 
by lattice invariance. 

The results in Fig. showing the serrated stress-strain curves obtained from 
the numerical simulations, evidence the bursty character of slip-mediated plastic 
deformation in the model. The mechanical response of the crystal occurs through 
stress-relieving dislocation avalanches at many scales. Hardening is absent when 
only one or few compatible slips are activated as in simulations I, II, III, where, 
after the initial elastic climb the stress reaches a yielding value, after which it 








Paolo Biscari et al. 



Fig. 3 (Color online) Microstructures developed in plastic flow. The thin black lines outline 
the undeformed reference configurations, (a) Band-like stress distribution during simulation I 
(section on the (a, m)-plane). (b) Detail of the stress peaks at the highest plastification zones of 
a shear band, with (c) the corresponding slip amplitudes, where double slip has also occurred, 
(d) Almost-coherent shear bands involving compatible shears in simulation II. Dark blue is 
unsheared; the other colors indicate the principal slips activated by the loading, and different 
combinations thereof, (e) Complex microstructure during simulation III, which involves mostly 
kinematically compatible slips, (f) Complex microstructure developed in simulation V, which 
activates mostly non-compatible slips. A similar result to V holds for simulation IV (not 
shown), which also involves mostly non-compatible slips. Dark blue is unsheared; the other 
colored domains involve various combinations of the fundamental slips, forming kinematically 
compatible as well as non-compatible pairs. 


remains substantially at a plateau value, which depends on the number and com¬ 
patibility properties of the involved slips. In contrast, hardening is present in the 
more general case of simulations IV and V (see the green and purple dashed lines 
in Fig.right), a typical effect in many experimental observations. As the initial 
crystal is fairly homogeneous, we observe in all cases the expected nucleation peak, 
common in other similar rate independent hysteretic phenomena, both theoreti¬ 
cally and experimentally The static defects have a marked influence 

on the first-cycle peak, which significantly decreases with increasing quenched dis¬ 
order (left inset in Fig. [^. The peak is anyway lowered and eventually eliminated 
in subsequent loading cycles due to the dislocation structures developed during 
the previous cycles (see the blue line in the left panel, and the red dotted line in 
the right panel of Fig. [^. 

Fig. [3] demonstrates the crucial role that compatibility among the involved 
principal slips plays in the developing microstructures, and in controlling the cre¬ 
ation of dislocations. Several features emerge, showing how the discrete plastic 
map captures many details of defect-mediated deformation. Panels (a)-(d), 
corresponding to loadings I and II, highlight the similarities in plastic processes 
which involve a single or several compatible slips. In both cases, at an intermediate 
yielding stage the plastic deformation is organized in bands which are in rough 
alignment with the optimal directions for stress relaxation at the discontinu- 
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iog,„(AE) iog^(,(AE) 

Fig. 4 Log-log plot of the heavy-tailed probability density of the energy drops along the 
plastic hysteresis (logarithmic binning). Left: data from simulations I (circles), II (squares), III 
(crosses). Right: data from simulations IV (diamonds) and V (circles). Data for simulations 
II, III, and V, are shifted down for clarity. The solid lines, drawn to guide the eye, indicate 
power laws with exponent 1.2 (left) and 1.5 (right). See also Table^ 


ities, with dislocations mostly concentrated at these non-fully-coherent interfaces. 
However, depending on the boundary conditions or the lattice orientation, much 
more complex microstructures may arise also in the case of activated slips which 
are mostly compatible, as happens in simulation III, see Fig.j^e). Despite the com¬ 
plexity of the microstructure developed in simulation III, the fact that it involves 
mostly compatible slips is likely at the origin of the lack of hardening observed in 
this case in the material response. Complex microstructures are a fortiori created 
in simulations IV and V, wherein the more generic loading conditions activate 
a (sub)set of non-compatible slips, see Fig. [^e) showing simulation V as an ex¬ 
ample. In these more general loadings the non-compatible fronts involve strong 
defect nucleation, although the deforming lattice typically attempts to accommo¬ 
date the boundary conditions through selective (cross-)shp combinations, gener¬ 
ating compatible interfaces as much as permitted by local minimization. However, 
the evolving structures are unable to fully obviate to the energetic cost imposed 
by non-compatibility, and the ensuing rapid accumulation of internal stresses and 
defects is a likely factor [8] in the hardening observed in the stress plateau of 
simulations IV and V, as in the green and purple dashed lines in Fig. i (right). 

Intermittency. Fig.j^shows the statistical properties, derived from the energy 
drops accompanying the slip avalanches during plastihcation, which generate the 
spatial patterns in Fig. The three loading conditions I, H, HI, in which only 
compatible shears are activated, produce largely equivalent heavy tailed statistics. 
We see that the size distributions associated to these fast dissipation events decay 
compatibly with power laws over about two decades, with exponents near or just 
below 1.2, see Table Simulation IV and V, on the other hand, produce heavy¬ 
tailed statistics consistent with a higher power-law exponent near 1.5. Analogous 
results hold for the stress drops, with respective exponents near 1.5 and 1.7, as 
is also reported in Table Interestingly, the higher exponents pertain to the 
cases showing evidence, such as hardening, of a higher defect density induced 
by the imposed loads due to the kinematic non-compatibility of the activated 
slips. This might lead to dislocation entanglements possibly screening long-range 
defect interactions and acting as barriers to free dislocation glide, with, in turn, a 
reduction of avalanche size and an increase of the statistics’ exponents, see also m- 
Detecting such compatibility-related effects, which are known to play important 
roles also in the (reversible) martensitic transformations of crystalline substances 
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Simulation 

Exponent energy drops 

Exponent stress drops 

I 

1.22 ±0.01 

1.46 ±0.03 

II 

1.19 ± 0.01 

1.49 ±0.04 

III 

1.15 ±0.02 

1.48 ±0.03 

IV 

1.54 ± 0.05 

1.75 ±0.05 

V 

1.57 ± 0.05 

1.72 ±0.04 


Table 1 Exponents, computed by the maximum likelihood method m, for the distributions 
of the energy and stress drops pertaining to simulations I-V. See also Fig.[^ 


such as shape-memory alloys [SHE], is a main advantage of our method based on 
an exact, 3D multiplicative evolution law § for the plastic strain. 

The above results on the dependence of the plasticity exponents on the specific 
wells involved in the flow, and thus, on the specific symmetry of the crystal and the 
specific loading imposed to the material, may be juxtaposed the analogous effects 
occurring in the martensitic transformations of crystalline materials such as shape- 
memory alloys. In the latter, it is known that scale-free behavior is characterized 
by exponents that are material- and loading-dependent, due to the symmetry, 
the compatibility properties, and the number, of the energy wells involved by the 
specific boundary conditions driving the transformation [iOlliTHiillSi] . Our present 
results evidence analogous effects also in scale-free crystalline plastic flow. 


4 Conclusions 

We have developed a nonlinear 3D phase field model accounting for the infinite 
and discrete invariance group of periodic lattices, which captures the main sig¬ 
nificant effects in crystal plasticity at microscopic scales. The non-compatible 
plastic distortion jumps resolve the single lattice defects, and plasticity in the 
model spontaneously proceeds through intermittent defect nucleation and evolu¬ 
tion in the lattice, with jump fronts orienting themselves to reduce the elastic 
energy, and with the plastic distortion always driven by the phase fields towards 
crystallographically-preferred configurations. Via the relation G = E~^GL(3,Z)E 
earlier recalled, and by adopting the appropriate set of principal slips Sa, the 
present approach can be adapted to investigate such phenomena also in crystalline 
materials whose symmetry is other than simple cubic as considered above. 

The numerical results show how both simple and complex strain topologies 
may develop in the deformed crystal, depending on the external loading and the 
kinematic compatibility properties of the different slip systems that are activated. 
Simulated plastic flow also shows a power-law decay for the avalanching events. The 
simplest cases mostly involving compatible active slips exhibit band-structured 
plastic deformation, low dislocation density and no hardening. For these we esti¬ 
mate low exponents agreeing with the low exponents recently determined in the 
detailed DDD simulations of plastic yielding [22 and in the re-analysis in |63l 
[26] of experimental microplasticity data in [^1191157] . We find higher exponents in 
more generic plastification conditions, when the subset of activated slips are mostly 
non-compatible, producing complex microstructures and intermittent yielding ac¬ 
companied by hardening. These results suggest that microscale crystal plasticity 
should not belong [29] to the universality class characterized by exponents com- 
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patible with mean-field depinning [6211161156] . In general, scale-free plastic flow in 
crystals may actually be material and loading dependent and lack universality, 
as is also advanced in recent work m, and in possible analogy to the material 
and loading dependent statistical features of avalanching behavior also observed 
in reversible martensitic transformation in crystals [iQ l ii H ii H iSi]. 
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contract. 
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